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Abstract. 

fj , We study a gas flow in the Laval nozzle, which is a convergent-divergent tube that 

O^' has a sonic point in its throat. We show how to obtain the appropriate form of the 

^ • tube, so that the acoustic perturbations of the gas flow in it satisfy any given wavelike 

^*J[ equation. With the help of the proposed method we find the Laval nozzle, which is an 

acoustic analogue of the massive scalar field in the background of the Schwarzschild 

black hole. This gives us a possibility to observe in a laboratory the quasinormal 

ringing of the massive scalar field, which, for special set of the parameters, can have 

O^ ' infinitely long-living oscillations in its spectrum. 

cn 
^' 

^ I 1. Introduction 

Massive fields in the vicinity of black holes have been studied during the last two decades 

'k/{ • (see [1] for review). It was found that their behaviour is qualitatively different from the 

H ■ behaviour of the massless fields. The response of a black hole upon any perturbations 

" " " at late times can be described by a characteristic spectrum of exponentially damped 

oscillations. The spectrum of a massive field, for some particular values of the 

parameters, has the oscillations with a very small decay rate in their characteristic 

spectra. These oscillations, that behave similarly to standing waves, were called 

quasiresonances [21 [3]. The asymptotical behaviour of massive fields is also different: 

one observes the oscillating tails that decay as inverse power of time, which is universal 

at the asymptotically late times [4]. 

Yet, since massive fields are short-ranged, we cannot expect the observation of 
their signal from black holes in near-future experiments. An attractive possibility 
for experimental study of the massive fields in the background of a black hole is a 
consideration of the acoustic analogue. This is a well-known Unruh analogue of a black 
hole [5], which is an inhomogeneous fluid system, where the perturbations (sound waves) 
can be described by a Klein-Gordon equation in the background of some effective curved 
metric [6]. 
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The sound waves in a fluid can propagate from a subsonic region to a supersonic 
one, but tliey cannot go back. Tlierefore, sonic points in a fluid with a space-dependent 
velocity form a one-way surface for the sound waves, which is called "acoustic horizon" 
by similarity with the event horizon of the black hole. 

The works of Unruh stimulated the study of various acoustic systems, such as: 

(i) the "draining bathtub" which is an analogue of a rotating black hole [3 El E]; 

(ii) the Bose- Einstein condensate p^ ITTl [T2l [13] . which, in the regime when the thermal 
fluctuations can be neglected, allows to observe a phonon analogue of the Hawking 
radiation [H [151 [Ml El [El [19] ; 

(iii) the so-called optical black holes [20] due to sound waves in a photon fluid of an 
optical cavity and inside an optical fiber [21]; and others [221 ESI [21]. 

Within the analogue-gravity approach one considers hydrodynamical equations as field 
equations in some effective background which is not a solution to Einstein equations. 
Using this approach the one- dimensional perturbations in the Laval nozzle were studied 
in [25]. It was found that the perturbations of the gas flow in the Laval nozzle can be 
described by a wave-like equation with the effective potential, which depends on the 
form of the tube. The inverse problem for the correspondence of the form of the Laval 
nozzle to the Schwarzschild black holes has been solved in [26], where the form of the 
Laval was found in order to obtain acoustic analogue for the perturbations of massless 
fields. 

Here we describe a method, which allows us to find an appropriate form of the 
Laval nozzle for any given effective potential. We use this method to obtain the 
acoustic analogue of the massive scalar field in the background of the Schwarzschild 
black hole. This paper is organized in the following form: In the section II we give the 
basic equations for a one- dimensional flow in the Laval nozzle and its perturbations. In 
the section III we describe the numerical method, which allows us to find the appropriate 
nozzle form in order to mimic any given effective potential. In the section IV we apply 
the method to find the form of the Laval nozzle, which is an acoustic analogue of the 
massive scalar field in the Schwarzschild background and show the corresponding time- 
domain profiles. Finally, in the conclusion, we discuss the obtained results and open 
questions. 

2. Basic equations 

A perfect fluid in the Laval nozzle can be described by the continuity equation and the 
Euler equation, that read, respectively, 

dtipA) + d,{pvA) = 0, (1) 

p{dt + v-V)v = - Vp, (2) 

where p is the density of a gas, v is the fluid velocity, p is the pressure, and A is the 
cross-section area of the nozzle. Following [25] we assume that the fluid is isentropic 
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and the pressure depends only on the density 

P oc p^, (3) 

where 7 is the heat capacity (7 = 1.4 for the air). 

Assuming that the flux is irrotational V x t; = 0, the velocity can be expressed 
as t; = V$, where $ = J vdx is the velocity potential, which satisfies the Bernoulli 
equation 

dt<l> + \ {d^^f + h{p) = 0. (4) 

We study linear perturbations of the flux, i.e. we consider the fluid density p and 
the velocity potential $ as 

p =-p + 5p, p'>\5p\ (5) 

$ = $ + 0, |9,<l|»|9,0| (6) 

where p, $ are the background dynamical quantities which satisfy ([1]) and dl]), 5p and 
describe the perturbations, which are considered small so that we neglect the higher- 
order corrections. 

We introduce the function H^{x), 

H^{x) = g^/^ I rfte^'^(*-'^(^)V(i, a;), (7) 

with 

pA , , f \v\dx ,„. 

g = —. <x) = / ^^, 8 

Cs -J C% — V^ 

where Cg is the sound speed, 

,,= /*=& (9) 

\ dp \j p 

We find that Hi^ satisfies the Schrodinger-type wave-like equation 

{^2+^"- F(x*)) H^x*) = (10) 






9 (Pg 1 ( dg ^ 
2 dx*"^ 4 I dx* 
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with respect to the new variable 

^^^rc^cjx (12) 

J cj — v^ 
where k = uj/csq and c^o is the stagnation sound speed. 

The coordinate x* is the tortoise coordinate for the analogue black hole: x* = —00 
at the throat and x* = 00 corresponds to the spatial infinity (x = 00). 

Following [26], we measure A and p, respectively, in the units of cross-sectional area 
at the throat {Ath) and the flux stagnation density (po), and choose the arbitrary factor 
for the function g in such a way that 

*=2^- A-^-il-p'^-^'Tp- (13) 
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Then the cross section area can be expressed as a function of g as 



V2 {2g^ (l - v/T^^ 



1/(7-1) 



A = ^ \ . (14) 



We find also that 



v^ 



,2 ^_i.2/(l-Vl-^-^)-l). (15) 

Since the gas velocity is equal to the sound velocity at the acoustic horizon, we obtain 

7 + 1 3 



9 



- r- . - ^- (16) 

horizon l^pi^^f^X VS 



3. The nozzle form from a given effective potential 

Linear perturbations of a spherically-symmetric black hole, after decoupling of the time 
and angular variables, can always be reduced to the following wave-like equation 
/ d 



^^^^^+u'-VinjRin = 0, (17) 

where the effective potential V = V{r*) depends on the parameters of the field and the 
black hole and the tortoise coordinate is defined as 

dr 



fir) 
where /(r) depends on the parameters of the black hole. 

In order to find the form of the Laval nozzle which is an acoustic analogue of the 
black hole perturbations we equate the tortoise coordinates and the effective potentials 
of the equations f lTO|) and f lT7|) 

f(r)f'ir)g'ir) + f{rfg'{rf _ f{rfg'{rf ^ 

2g{r) Ag{ry ^''^^ ^ ^ 

From dx* = dr* and equations flT^ and flT^ we find the relation between the coordinate 
of the nozzle and the radial coordinate of the metric r: 



dx = (^+^-'^^;^'(^-V^-^(^H)^^ (20) 

/(r)(7 - l)^2g{ry{l-^l-g{r)-^) 

li g{r) is known, from the equations ( IT4|) and ( l20l) . one can find the function A{x), which 
describes the nozzle form. 

In order to find g{r) we make the substitution g{r) = h{rY. Then the differential 
equation (IT^ reads 

f{rfh'\r) + f{r)f\r)h\r) - V{r)h{r) = 0. (21) 

Since the function /(r) vanishes at the event horizon r = r+, the linear equation 
(1^ always has a regular singular point there. Using the Frobenius method we expand 
the general solution to the differential equation near the event horizon as 

/i(r) = Ci/ii(r) + C2/i2(r), (22) 
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where Ci and C2 are arbitrary constants, 

h,{r) = (r - r+)^i ( 1 + £ a^^ir - r+)A , (23) 



n=l 



hir) = hi{r) ln(r - r+) + (r - r+)^2 ^ 6„(r - r+)", 

ra=0 

when Ai — A2 is an integer, and 

h,{r) = (r - r+)^^ U + £ &n(r - r+)M , 

otherwise, A2 < Ai are the roots of the indicial equation and depend on the given 
functions /(r) and V{r). 

In order to satisfy (1161) . one of the roots must be zero. f'{r^) > iniphes that the 
other root is negative. Hence, for A2 < Ai = 0, h2{r) is always divergent at the horizon 
r = r^ and we choose C2 = 0. Therefore, from flT^ we find that 



7 + 1 
Cl = 



We expand ( 123|) near the event horizon and find h'{r+), which completely fixes the 
initial value problem at r = r+. Then, we are able to solve numerically the equation 
fl2T]) using the Runge-Kutta method for r > r+. 

4. Acoustic analogue for the massive scalar field 

We consider the massive scalar field in the background of the Schwarzschild black hole, 
given by the line element 

Hr'^ 2M 

ds^ = f(r)dt^ - -77-T - r\de^ + sin^ ed(P^), fir) = 1-^1^, (24) 

fir) r 

where M is the mass of the black hole. Hereafter we measure all the quantities in units 
of the black hole horizon, i.e. r^ = 2M = 1. 

The scalar field \1/ satisfies the Klein-Gordon equation 

(V'V^ + m^)m = 0, (25) 

where V^ is the covariant derivative, m is the field mass. The equation (l25l) in the 
background (121|) reads 
1 



d^. ( g^''^\g\d^'^ ) + m'^ = 0. (26) 

After the separation of the angular and time variables 

^{t,r,e,<p) = ^Yr{e,<P)e-'''' (27) 
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we obtain the wave-like equation (1171) with the effective potential 

nr) = /(r)(^ + ^+m^). (28) 

In order to show the time-domain evolution of perturbations we use the 
discretization scheme proposed by Gundlach, Price, and PuUin [27]. We consider the 
time- dependent equation 

( cP d^ \ 

- — -\/(r)U(t,r*) = 0. (29) 



\dr*'^ dt"^ 

Rewriting f l29|) in terms of the light-cone coordinates du = dt—dr* and dv = dt+dr*, 
we find that 

$(Ar) = $(H^) + $(E) - $(S) - — V'(^) [^{W) + $(E)] + C(/i^), (30) 

8 

where the point N, M, E and S are the points of one square in a grid with step h in the 

u-v plane, as follows: S = {u,v), W = {u + h,v), E = {u,v + h) and N = {u + h,v + h). 

With the initial data specified on two null-surfaces u = uq and f = f o we are able to 

find values of the function \l/ at each of the points of the grid. 

On the figures ([1]) and ([2]) we show the forms of the nozzle that are acoustic analogs 
of the massive field with particular values of the mass for which nearly infinitely long- 
living oscillations exist. These oscillations called quasiresonances one can observe on 
the corresponding time-domain profiles. We see that in the tube of a particular form 
the decay rate of sound waves is almost zero for some tone which is an analogue of the 
quasiresonance of the massive scalar field. 

Although we present here only the nozzles where the quasiresonances can be 
observed, the method described above can be used to construct an analogue for any 
finite mass of the in such a way that the sound waves in the nozzle will have the same 
behaviour as the massive scalar field in the background of the Schwarzschild black hole. 
From the figures ([T]) and ^ one can observe that the higher mass is, the quicker the 
nozzle cross-section grows, diverging at the end. However, as it was pointed out in [26] , 
this does not lead to a problem with the presented model because of the freedom of the 
choice of the units of length. One can rescale the nozzle along the transversal axis in 
order to make the cross-section change as slowly as one wants. This change of the scale 
changes proportionally the frequencies of the sound in the nozzle. 



5. Conclusion 

We have considered the Laval nozzle as an acoustic analogue of the massive scalar field 
in the background of the Schwarzschild black hole. We presented the general method to 
determine the form of the Laval nozzle, such that the sound waves in it are described 
by a given effective potential, what can be used to study an analogue of black hole 
perturbations in a laboratory. The method can be used to obtain the forms of the 
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Figure 1. The form of the Laval nozzle for the spherically symmetric {£ — 0) massive 
scalar field (a): m = 0.8 (green, narrow), m — 2.84 (blue), and m = 4.86 (magenta, 
wide), and the corresponding time-domain profiles: (b) m — 0.8, (c) m — 2.86, (d) 
TO = 4.86. 




Figure 2. The form of the Laval nozzle for the £ ^ 1 massive scalar field (a): 
TO. — 1.06 (green, narrow), to. — 2.30 (blue), and m — 4.40 (magenta, wide), and 
the corresponding time-domain profiles: (b) to = 1.06, (c) to = 2.30, (d) to = 4.40. 
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nozzles, which are acoustic analogs of other spherically symmetric black holes. The 
acoustic analogs for perturbations of Reissner-Nordstrom(-de Sitter) black holes and 
their higher- dimensional generalizations, black strings, and Gauss-Bonnet black holes 
are of special interest. For some set of the parameters the higher- dimensional black 
holes and black strings suffer instability [2H1 EHl [30] , which in the corresponding nozzle 
can manifest itself as increasing of the sound amplitude. It is clear that for a large 
amplitude of the sound waves cannot be described within the linear approximation 
and the considered analogue between the linear perturbations cannot be apphed. 
Nevertheless, we believe that the consideration of different physical systems, which have 
linear instability in the same parametric region, could help us to understand better its 
nature. 
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